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Abstract 

We present a density functional scheme for calculating the pair density (PD) by means of the 
correlated wave function. This scheme is free from both of problems related to PD functional 
theory, i.e., (a) the need to constrain the variational principle to iV-representable PDs and (b) 
the development of a kinetic energy functional. By using the correlated wave function, the 
searching region for the ground-state PD is substantially extended as compared with our previous 
theory [Physica B 372 (2007), in press]. The variational principle results in the simultaneous 
equations that yield the best PD beyond the previous theory, not to mention the Hartree-Fock 
approximation . 



Note: The previous paper [Physica B 372 (2007), in press] can be downloaded from "Arti- 
cle in press" in the website of Physica B. 
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I. INTRODUCTION 



The pair density (PD) functional theory has recently attracted particular interest because 



it provides the obvious way to improve on the density functional theory (DFT 



it provides tJ 

w. Zi 



]j_and then 



Ziesche first proposed the PD functional theory about a decade ago 
many workers followed his work and have developed a variety of approaches 

Very recently, we have proposed an approximate scheme for calculating the PD on the 
basis of the extended constrained-search theory [n], [ll|, Q Q • By introducing the 



noninteracting reference system 







11 1, the resultant PD corresponds to the best solution 



within the set of the PDs that are constructed from the single Slater determinant (SSD). 
This PD functional theory has two kinds of merits. The first one is that the reproduced PD 
is necessarily iV-representable. This is a strong merit because the necessary and sufficient 
conditions for the iV-representable PD have not yet been known 2, 
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The second merit is the tractable form of the kinetic energy functional. 
The kinetic energy functional cannot exactly be written by using the PD alone. Some 
approximation is required to be introduced [tJ. In this theory, we have successfully given an 
approximate form of the kinetic energy functional with the aid of the coordinate scaling of 

electrons MM- 

n 

On the other hand, we also have the remaining problem in it|10j. Namely, there exists the 
possibility that the solution might be far from the correct value of the ground-state PD. This 
is because the searching region of the PDs may be smaller than the set of iV-representable 
PDs. In order to improve the PD functional theory with keeping the above-mentioned 
merits, we have to extend the searching region of the PDs to the set of TV-represent able PDs 
as closely as possible. At least, we had better extend the searching region beyond the set of 
the SSD- represent able PDs 27]. 

In this paper, we shall employ the strategy for reproducing the PDs not by means of 
the SSD, but through the correlated wave function. As the correlated wave function, we 
adopt the Jastrow wave function that is defined as the SSD multiplied by the correlation 



function 



28 



291 ]. Owing to the correlation function, the searching region substantially be- 



comes larger than the set of the SSD-representable PDs. Of course, the reproduced PD is 
kept iV-representable because the PD is calculated via the Jastrow wave function that is a 
kind of antisymmetric wave functions. Also the second merit is not missed in the present 
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scheme by utilizing the result of the scaling property of the kinetic energy functional. 

Organization of this paper is as follows. In Sec. II, we provide the preliminary definitions 
of various quantities that appear in the following sections. In Sec. Ill, by means of the 
variational principle with respect to the PD, we derive simultaneous equations that yield 
the best PD that is superior to the previous oneficl]. Such equations are quite tractable, the 
computational method of which are also proposed in Sec. III. 



II. PRELIMINARY DEFINITIONS IN THE PD FUNCTIONAL THEORY 

In this section we give the preliminary definitions which will be used in the present 
scheme. The PD is defined as the diagonal elements of the spinless second-order reduced 
density matrix, i.e., 



7^(rr'; rr') = (\&| - // $ + (r, ', 7/)$(r', r]')ijj(r, r])dr]dr]' \^f) 
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(1) 



where tp(r, rj) and "0 + (r, rj) are field operators of electrons, and \l/ is the antisymmetric wave 
function, and r and rj are spatial and spin coordinates, respectively. We shall consider a 
system, the Hamiltonian of which is given by 



H = T+W+ 



j drp(r)v ext (r), 



(2) 



where T, W and p(r) are operators of the kinetic energy, electron-electron interaction and 
electron density, respectively, and v ex t(r) stands for the external potential. In the similar 
way to the extended constrained-search theory }l2. 13 . 1^ . the universal functional is defined 

as 



Min (V\T + W\V) 

Vlr—Kyp) ( rr ';rr') 



T + W 



vf[7 {2) ]), 



(3) 



where ^ — > 7^ 2 - ) (rr'; rr') denotes the searching over all antisymmetric wave functions that 
yield a prescribed 7 < - 2 - ) (rr / ; rr'). In the second line, the minimizing wave function is expressed 
as ^/[7^ 2 ^]. By using Eq. (j3J), the Hohenberg-Kohn theorems for the PD functional theory 
can be easily proved|l|, Here we show only their results jlo|: 



(4) 
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and 

E = Min£[ 7 {2) ] 

*y(2) 

= £[7$ 2) ], (5) 



(2) 

where fy 0j E and % are the ground-state wave function, ground-state energy and ground- 
state PD, respectively, and where E^ 2 -*] is the total energy functional that is given by 

£[ 7 (2)] = F[y™] + -3-r \\ drdrV^(r)7 (2) (rr';rrO. (6) 



N-l 

Equations (@J) and (jSJ) correspond to the first and second Hohenberg-Kohn theorems, respec- 
tively. Let us suppose that 

T[^] = (y[ 1 W}\f |*[ 7 ( 2 )]), (7) 

then Eq. (JUJ) is rewritten as 

£[7 (2) ] = T[ 7 (2) ] + e 2 // drdr'^^^ + J J drdr'^(r) 7 ( 2 )(rr'; rr'), (8) 

where, in the second term, we use the fact that the expectation value of W is exactly 
written in terms of 7 ^ 2 ^(rr'; rr'). Equation (jHJ) is the starting expression for the total energy 
functional in the PD functional theory. 

As mentioned in Sec. I, the kinetic energy of the PD functional theory cannot be exactly 
expressed by the PDs alone. In other words, we have to employ the approximate form in Eq. 
(JHJ). So far, the kinetic energy functional of the PD functional theory has been developed 



by several workers 
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30l |3l| . In this paper, we make use of an approximate form of the 
kinetic energy functional which has been derived by utilizing the scaling property of the 
kinetic energy functional [lfl Isil]. The explicit form is given by 

T[ 7 ( 2 )] = K jj drdr' 7 (2) (rr'; rr') s , (9) 

where K is an arbitrary constant. 

III. SINGLE-PARTICLE EQUATIONS 

Equation (JSJ) corresponds to the variational principle with respect to the PD. The search- 
ing region of the PDs should be of course within the set of iV-representable PDs. For that 
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purpose, we shall introduce the searching region of the PDs that are calculated from the 
correlated wave functions. The searching region is substantially extended as compared with 
the previous theory [l^, because it is restricted within the set of SSD-representable PDs. 
Extension of the searching region can be regarded as one of appropriate developments of the 



PD functional theory 



27]. 



In this paper we adopt the Jastrow wave function as the correlated wave function. The 
explicit evaluation of the PD using the Jastrow wave function is actually very hardj3, S|- 
As a consequence, the approximation technique to evaluate the PD has been developed espe- 
cially in the field of nuclear physics. The expectation value of the PD operator with respect 
to the Jastrow wave function can be systematically expressed with the aid of the Yvon- 



Mayer diagrammatic technique 



28, 



291 ] . Here we shall use the lowest-order approximation of 



the expectation value of the PD operator. 
The Jastrow wave function is defined as 



28, 



29] 



%j(xi,x 2 , ,x N ) = — = J[ f(r i:j )® S sD(xi,X2, ,x N ), (10) 

V^N l<i<j<N 

where &ssd(%i, %2, ,%n) is the SSD, and where /(r^) = / (|r-j — r,]) is the correlation 
function, and where is the normalization constant. Suppose that the correlation function 
is chosen to satisfy the cusp condition for the antisymmetric wave function. The lowest-order 
approximation for the expectation value of the PD operator is given by |28j 

7 ( 2 )(rr';rr') = \f (|r - r'|)| 2 7 § D (rr'; rr'), (11) 

where ^^(rr'; rr') is the expectation value of the PD operator with respect to the SSD. 

Supposing N orthonormal spin orbitals of the SSD are denoted as {i/j^(x)}, then Eq. (fTTj) 

is explicitly expressed as 

1 N f f 

7 (2) (rr / ;rr) = - ]/ (|r - r'|)| 2 £ // drjdr]' {^(^^(x')^^)^^') 



£11,^2 = 1 



Next, let us consider the variational principle with respect to the PD, i.e. Eq. (jSJ). The 
variation of the PD is performed via the spin orbitals of Eq. ()12|) with the restriction that 
they are orthonormal to each other. Using the Lagrange method of undetermined multipliers, 
we minimize the following functional without the restriction: 

= £[ 7 (2) ] -E^{/ V£(z)^(x)dz-<v}, (13) 



where Eqs. (JHJ), © an d (|12|) are used in the first term on the right-hand side. The minimizing 
condition 5fl [{VVt}] = immediately leads to 

/ ^{^{xij^x^^x) - 



x 1/ (|r - ri |)| 2 |^7 (2) (rr i; rn)* + + j^-j (<w(r) + t^n))} 

= X^mW^Oz), (14) 

V 

where the chain rule for the functional derivatives is utilized. The Lagrange multipliers e^ v 
should be determined by requiring that the spin orbitals are orthonormal to each other: 



ip* fl (x)ijj l/ (x)dx = 8^. (15) 

In the similar way to the previous theory [lo|. we can simplify the above equations by means 
of a unitary transformation of the spin orbitals. It is easily shown that e M „ forms the 
Hermitian matrix. Suppose that the unitary matrix which diagonalizes is written by 
U^, then 

is satisfied, where is the diagonal element of the diagonal matrix. Let us consider the 
following transformation of the spin orbitals: 

VVi(aO = J2 U ^Xv(%)- (17) 
Substituting Eq. (fTTj) into Eq. (JHJ), and using Eq. (fTTIjl . we obtain 

v J 

x 1/ (|r - ri |)| 2 |^ 7 (2) (rri; m)* + + ^ (<w(r) + <w( ri ))} 

= e^Xiiix). (18) 
Also, Eq. (|15jl is transformed into 

/ Xl t (x)xu(x)dx = 8^. (19) 

Here note that the expression for 7^) (it'; it') in Eq. (|18|) is kept invariant under the unitary 
transformation. This is confirmed by substituting Eq. (|17|) into Eq. (|12|). i.e., 

1 N r r 

-|/(|r-r'|)| 2 *Vdv' {x;MK 2 (x')x»Ax)x» 2 (x') 

- xlMxUx')x^)x^')} ■ (20) 



7^ (it'; rr' 



Ml>A»2=l ' 



Equations (fTSj) and (|P3J) are the simultaneous equations, and the solutions yield the best 
PD within the set of PDs that are calculated from Eq. (|20|) . 

Our previous work may be the first proposal of a computational approach that deals with 
problems related to the PD functional theorv[lo|. The present scheme is also a computational 
approach, and further improves on the previous theory concerning the searching region of 
the PDs. In that sense, it would be useful to consider a computational procedure for solving 
the simultaneous equations (fTHj) and (fTUj) . 

The procedure proposed here is similar to that of the Hartree-Fock equation j^. In order 
to make the computational procedure readily comprehensible, let us rewrite Eq. (fTHjl as 



{F(r)-i 5 } X 6(x) = Gs(x) (21) 



with 



N 

F(r)= / dx 1 \f(\v-v 1 \)\ 2 J2\Xu(xi 

J v=l 



x j^7 (2) (rri;rr 1 )U^^ + -^-(^(r) + ^(r 1 ))), (22) 
J r — rJ I\ — 1 I 



f 2 { N 

G s (x) = I dx 1 \f (|r - r'|)| 2 I ^ xt{ x i)xs{xi)Xu 



x 



x j^7 (2) (rri;r ri )i + ^ f + ^—r (v ext {r) + Uext(ri))) , (23) 

{ 3 |r — ri| N — 1 J 

where the spin orbital Xs{ x ) is the solution of Eq. (J21j) . and should be determined in a 

self-consistent way. Here note that the right-hand side of Eq. (}2*Tj) comes from the second 

term of Eq. (|20|1. and explicitly depends on the spin orbital Xs( x )- The key point to get the 

self-consistent solution is that spin orbitals of the previous iteration are used in calculating 

F(r) and G$(x) 32}. By solving simultaneously Eqs. (fT^j) and (j2~Tj) with this technique, we 

can get a new set of spin orbitals and energy parameters sVs. We continue such a procedure 

until the self-consistency for the solutions is accomplished j^. 

IV. CONCLUDING REMARKS 

In this paper, we propose the PD functional theory that yields the best PD within the set 
of PDs that are constructed from the correlated wave functions. Compared to the previous 



one 
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11 1, the present theory has the following features. 



1. The present theory is superior to the previous one in that the searching region of 
the PDs is certainly larger than the set of SSD-representable PDs without missing the 
merits of the previous theory Q]. This means that the resultant PD is more reasonable 
than that of the previous theory. 

2. The predominance of the present scheme can also be shown from the viewpoint of the 
total energy. If the correlation function is chosen to be unit, then the present theory is 
reduced to the previous one exactly. It has been already proved that the total energy 



of the previous theory is better than that of the Hartree-Fock approximation [11 1. 
If the correlation function is chosen most appropriately, then the searching region 
is substantially equivalent to the set of PDs which are calculated by varying both 
correlation function and spin orbitals in Eq. (|20|). Therefore, the total energy of the 
present scheme is necessarily more sound than the previous onejlO^. and needless to 
say, than that of the Hartree-Fock approximation. 

3. In addition to the above merits, the present scheme has the feature that deserves 
special emphasis. Due to the fact that the PD functional theory is still a developing 
field, there hardly exist the computational approaches so far. Our previous work is 
perhaps the first paper to propose a computational approach that incorporates both 
of problems related to PD functional theory jiol. fll | . The present scheme is also a 
computational approach. The resultant simultaneous equations are quite tractable, 



111 ]. Also from such a viewpoint, the present scheme 



as well as the previous one[l(l 
seems to be valuable. 

Thus, the resultant simultaneous equations (|18|) and (I19|) yield the PD which is definitely 
closer to the ground-state PD than the previous theory |]J| . Next step is to perform the actual 
calculations so as to confirm to what extent the present scheme covers the iV-representable 
PDs. 

Finally, we would like to comment on the future prospect of the present theory. Although 
the present scheme utilizes the lowest-order approximation of the expectation value of the 
PD operator, the higher-order corrections can proceed systematically with the aid of the 
Yvon-Mayer diagrams 



28l Q. Of course, it is anticipated that the equations will become 
more complicated. But, from the methodological point of view, it is important that the 
theoretical framework has the potentiality to improve the approximation systematically. 
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